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Abstract: Gross primary productivity (GPP) of vegetation is an important constituent of the terrestrial 
carbon sinks and is significantly influenced by drought. Understanding the impact of droughts on 
different types of vegetation GPP provides insight into the spatiotemporal variation of terrestrial carbon 
sinks, aiding efforts to mitigate the detrimental effects of climate change. In this study, we utilized the 
precipitation and temperature data from the Climatic Research Unit, the standardized precipitation 
evapottanspiration index (SPET), the standardized precipitation index (SPI), and the simulated vegetation 
GPP using the eddy covariance-light use efficiency (EC-LUE) model to analyze the spatiotemporal change 
of GPP and its response to different drought indices in the Mongolian Plateau during 1982-2018. The 
main findings indicated that vegetation GPP decreased in 50.53% of the plateau, mainly in its northern 
and northeastern parts, while it increased in the remaining 49.47% area. Specifically, meadow steppe 
(78.92%) and deciduous forest (79.46%) witnessed a significant decrease in vegetation GPP, while alpine 
steppe (75.08%), cropland (76.27%), and sandy vegetation (87.88%) recovered well. Warming aridification 
ateas accounted for 71.39% of the affected areas, while 28.53% of the areas underwent severe 
aridification, mainly located in the south and central regions. Notably, the warming aridification areas of 
desert steppe (92.68%) and sandy vegetation (90.24%) were significant. Climate warming was found to 
amplify the sensitivity of coniferous forest, deciduous forest, meadow steppe, and alpine steppe GPP to 
drought. Additionally, the drought sensitivity of vegetation GPP in the Mongolian Plateau gradually 
decreased as altitude increased. The cumulative effect of drought on vegetation GPP persisted for 
3.00-8.00 months. The findings of this study will improve the understanding of how drought influences 
vegetation in arid and semi-arid areas. 
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1 Introduction 


With the escalating global warming crisis, mitigation measures such as carbon emission reduction 
and the promotion of carbon neutrality have attracted worldwide attention, leading to 
international consensus on the adoption of green and low-carbon strategies (Ji and Lin, 2022; 
Zhao et al., 2022). Central to these discussions is the role of vegetation gross primary productivity 
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(GPP), which quantifies the amount of organic carbon content captured through photosynthesis in 
a unit of time. It determines the amount of initial materials and energy that enter the ecosystem 
(He et al., 2022b). As an important component of the global carbon sink, vegetation GPP 
profoundly affects the terrestrial carbon cycle. 

Long-term and high-precision GPP series are prerequisites for studying the spatiotemporal 
dynamics and influencing factors of vegetation productivity. Currently, vegetation GPP is mainly 
estimated through ground surveys, flux observation, remote sensing inversion, and model 
simulation. Traditional GPP monitoring methods, such as ground survey and flux observation, are 
primarily based on a sample scale. Due to the strong spatial heterogeneity of GPP, direct GPP 
observation based on a sample scale cannot reflect the changes in large-scale continuous 
vegetation productivity; consequently, conducting a comprehensive analysis of the spatiotemporal 
differentiation and variability characteristics of GPP at a regional scale would be challenging 
(Piao et al., 2013, 2022). In recent years, remote sensing inversion and model simulation methods 
have developed rapidly due to their capacity to characterize large-scale vegetation productivity 
effectively. Hence, they are widely used to calculate regional-scale vegetation productivity (Xiao 
et al., 2019; Xie and Li, 2020; Chu et al., 2021; Pan et al., 2021). For example, Zhang et al. 
(2020) used moderate-resolution imaging spectroradiometer (MODIS) GPP to analyze the 
spatiotemporal distribution and variation characteristics of vegetation GPP in China from 2003 to 
2018. The models commonly used for GPP estimation include machine learning models (based on 
flux tower data), ecosystem process models, and light use efficiency (LUE) models (Jiang and 
Ryu, 2016; Zhen et al., 2018; Wang et al., 2021b; Guan et al., 2022). Of these, LUE models 
require only a few input parameters. Moreover, they can directly obtain some key parameters 
through the inversion of remotely sensed data. These models mainly employ the vegetation index 
to estimate the photosynthetically active radiation absorbed by the canopy and calculate light 
energy utilization efficiency, which is the optimal measure of photosynthesis in the vegetation 
canopy. Therefore, it is widely considered the most effective method to monitor large-scale 
vegetation productivity and is widely used to study spatiotemporal variation (Vanikiotis et al., 
2021; Pei et al., 2022; Xie et al., 2023). For example, Bo et al. (2022) analyzed the spatiotemporal 
variability and trends concerning GPP in China from 1982 to 2015, primarily relying on the 
optimum LUE (LUEopt) product. 

Additionally, some scholars have further enhanced LUE models by introducing additional 
factors to improve its accuracy. For instance, Zheng et al. (2020) improved long-term global GPP 
estimation by adding environmental variables. Yuan et al. (2010) developed the eddy 
covariance-light use efficiency (EC-LUE) model by substituting the Bowen ratio with the 
clearness index. This modification not only improved the accuracy of LUE models but also 
succeeded in developing a series of long-term GPP products. The EC-LUE model can be applied 
across various types of ecosystems and geographical regions because of its ability to map GPP 
over extensive areas with high accuracy and applicability while avoiding regional model 
parameterization (Yuan et al., 2010; Cui et al., 2017; Yu et al., 2018; Wang et al., 2020). 

Recent research has highlighted climate change as the primary contributor to vegetation GPP 
trends in various regions. Wang et al. (2022a) pointed out that changes in water conditions are 
among the dominant factors determining vegetation carbon sinks. Over the past five decades, 
global warming in general, coupled with the continuous decrease in precipitation, has led to a 
marked increase, globally, in the area covered by arid land (Hang et al., 2012; Feng and Fu, 
2013), and extreme drought events have become more frequent, bringing unprecedented 
challenges to the ecosystem. As a crucial factor affecting GPP, drought has gained the attention of 
scholars (Chen et al., 2021; Xiong et al., 2022). Drought, a natural phenomenon caused by 
prolonged water scarcity, affects the ecological carbon cycle by altering the underlying 
physiological processes, such as photosynthesis and respiration (Zhang et al., 2020). 

The commonly used drought evaluation indices include the Keetch-Byram drought index 
(Charlton et al., 2022), standardized precipitation index (SPI) (Hayes et al., 2011), self-calibrating 
Palmer drought severity index (Huang et al., 2022), and standardized precipitation 
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evapotranspiration index (SPEI) (Vicente-Serrano et al., 2010a). Focusing on the semi-arid 
Caucasus Ecoregion, Akhalkatsi (2017) analyzed the factors affecting crop growth and found that 
seed quality diminishes under soil water stress in arid years, and prolonged drought can lead to 
crop wilting and death. Wang et al. (2022b) used four different GPP datasets and SPEI data to 
quantify GPP reductions in global terrestrial ecosystems caused by varying drought severities. 
The study found that drought-induced GPP reductions intensify with the severity gradient. 
Additionally, grassland and cropland were found to be the ecosystems most vulnerable to 
increasing drought severity. Regarding the duration of adverse conditions, even a few months of 
drought is likely to adversely affect vegetation productivity. Precipitation deficits over a period as 
short as three months or even shorter could lead to severe drought (Wang et al., 2016). Peng et al. 
(2019) investigated the cumulative and lagged effects of drought across extra-tropical ecosystems 
in the northern hemisphere and found 1.00—4.00 accumulated months as a significant duration for 
the cumulative effect and 2.00—6.00 lagged months for the lagged effect. Many studies have 
emphasized the importance of considering the cumulative effect duration (CED) when assessing 
the impact of drought on vegetation GPP, especially in arid and semi-arid regions. For example, 
Wei et al. (2022) conducted a comparative analysis of the accumulated and lagged effects of 
drought on global grassland productivity and found the former to be stronger than the latter in 
most global grasslands. Zhan et al. (2022) employed the multi-scale drought and vegetation index 
datasets to analyze the impact of drought on vegetation in Yellow River Basin of China, 
concluding that vegetation areas with low water availability are more susceptible to the 
cumulative impact of drought. In addition to reduced precipitation, temperature rises significantly 
contribute to drought. Notably, "hot drought", which occurs when high air temperatures coincide 
with precipitation deficits, threatens ecosystems as it increases the stress impact caused by greater 
vapor pressure scarcity. Nevertheless, the extent to which hot and dry atmospheric conditions 
amplify the impact of inadequate precipitation on the ecological carbon cycle remains unclear 
(Dannenberg et al., 2022). 

The Mongolian Plateau is located in a transitional zone between arid and semi-arid climate 
zones and is home to various vegetation types that play an important role in the global carbon 
cycle (Li et al., 2023). In recent decades, warmer and drier climates have caused obvious soil 
drying and increased degradation of vegetation productivity, threatening the stability of regional 
ecosystems and the achievements of vegetation construction (He et al., 2022a). However, recent 
studies concerning the plateau's vegetation GPP have grown in number and achieved valuable 
results (Chen et al., 2022). Historically, scholars primarily used normalized difference vegetation 
index (NDVI) to evaluate the vegetation dynamics in this region, focusing on its response to 
climatic factors such as precipitation and temperature (Guo et al., 2021; Kang et al., 2021; Na et 
al., 2021). Gu et al. (2022) noted that since 1998, the cumulative effect of drought (7.00—12.00 
months) has been the main contributing factor to NDVI as opposed to the lagged effect 
(10.00-12.00 months). Yin et al. (2023) pointed out that moisture management must be addressed 
as the main challenge in ecological restoration processes in the Mongolian Plateau. Recently, 
Meng et al. (2023) identified water availability, temperature, and solar radiation as the most 
crucial climatic factors contributing to GPP change in 51.21%, 15.51%, and 18.63% of the 
Mongolian Plateau, respectively. Therefore, it is necessary to analyze the response of vegetation 
GPP to drought in the Mongolian Plateau. 

Against this backdrop, this study aimed to explore the spatiotemporal variability of GPP and its 
response to different drought indices in the Mongolian Plateau from 1982 to 2018. The impacts of 
a warming climate on the interaction between vegetation GPP and drought were also revealed. To 
achieve these objectives, we first identified the spatiotemporal characteristics of changing GPP 
trends and examined their sustainability using the EC-LUE-generated GPP data. Furthermore, the 
drought indices (SPEI and SPI) were utilized to analyze the drought trends concerning the 
Mongolian Plateau over the past 37 years. Finally, the effects of different drought indices on 
vegetation GPP were determined, and different effects of droughts due to climate warming were 
also discussed. The novelty of this study lies in the subdivision of the study areas into warming 
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aridification zones and non-warming aridification zones by analyzing the change in dry and wet 
trends of the two drought indices (SPEI and SPI) and quantifying the possible effects of 
temperature rise on the response of GPP to drought. This study provides a scientific basis for 
comprehensively evaluating the risk of GPP reduction under the single and combined influences 
of heat stress and drought stress. 


2 Study area and data 


2.1 Study area 


The "Mongolian Plateau" generally refers to the plateau region of Northeast Asia (37°22'— 
53°23'N, 87°43’-126°04'E), encompassing the entire territories of Mongolia and the Inner 
Mongolia Autonomous Region of China (hereafter referred to as "Inner Mongolia") (Luo et al., 
2021). It extends to the Da Hinggan Ling Mountains in the east, the Altay Mountains in the west, 
the Sajan and Kente Mountains in the north, and the Yinshan Mountains in the south (Yin et al., 
2022). Geographically, the Mongolian Plateau's terrain is elevated in the west and low-lying in the 
east (Fig. la), with an average elevation of 1580 m. The plateau straddles the temperate 
continental climate zone, receiving an average annual precipitation of about 200 mm; however, it 
can be higher than 400 mm in the northern and eastern mountainous areas and less than 100 mm 
in the southwest Gobi desert, as the distance from the ocean increases. Additionally, the 
temperature in the humid areas of the plateau is lower, while the temperature in the arid areas is 
higher. The lowest temperature in winter can reach —-40°C, and the highest temperature in summer 
can surpass 35°C (Zhang et al., 2023). 

The Mongolian Plateau hosts a complex and diverse range of vegetation. Except in the Gobi 
desert in the southwest, coniferous forest, deciduous forest, meadow steppe, typical steppe, desert 
steppe, alpine steppe, shrub, and sandy vegetation are distributed alternately across the plateau 
(Fig. 1b). Natural grassland, an important component of the temperate grasslands of the northern 
hemisphere, is the dominant vegetation type in this region. Due to climate change and the 
resultant decrease in water supply, the Mongolian Plateau often faces droughts of varying severity 
(Zhang et al., 2022b). 
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Fig.1 Digital elevation model (DEM; a) and vegetation type (b) of the Mongolian Plateau 


2.2 Data 


This study included GPP data from 1982 to 2018, sourced from the National Earth System Science 
Data Center (http://www.geodata.cn). This data were simulated by the EC-LUE model, with a 
temporal resolution of 8 d and a spatial resolution of 0.05° (Yuan et al., 2010). The model relied on 
only four variables: NDVI, photosynthetic active radiation, air temperature, and the Bowen ratio of 
sensitive to latent heat flux. For this product inspection, the average GPP product of 8 d across all 
years at 95 sites was extracted. When compared with the measured GPP in the flux sites, the 
maximum average relative error of GPP was 11.8%. The Global Land Surface Satellite (GLASS) 
products, characterized by long time series, high resolution, and high precision, offer a reliable 
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basis for the study of global environmental change. They can be widely used in the global, 
intercontinental, and regional atmospheres, vegetation cover, and water bodies of dynamic 
monitoring, and they can be combined with climate change characterization parameters such as 
temperature and precipitation and applied to the global change analysis. The precipitation and 
temperature details were obtained from a monthly time series dataset provided by the Climatic 
Research Unit (CRU) (https://crudata.uea.ac.uk/cru/data/hrg/) of the University of East Anglia, 
with a spatial resolution of 0.50°. SPEI data were obtained from SPEI base v.2.7 dataset. It was 
calculated based on the precipitation and potential evapotranspiration data provided by CRU TS 
4.05, with a time scale of 1.00-12.00 months and a spatial resolution of 0.50°. The precipitation 
data from CRU TS 4.05 were further used to derive SPI-related findings. The time scale of SPI is 
consistent with that of SPEI. The digital elevation model (DEM) data were sourced from the 
Geospatial Data Cloud (http://www.gscloud.cn). We processed the spatial resolution of the data to 
0.50°. The vegetation type map was obtained through the vectorization of 1:1,000,000 vegetation 
type maps of Inner Mongolia and Mongolia, allowing us to divide the vegetation types in the 
Mongolian Plateau into nine categories: coniferous forest (7.24%), meadow steppe (8.98%), shrub 
(2.33%), typical steppe (32.50%), deciduous forest (2.64%), alpine steppe (3.04%), cropland 
(4.66%), desert steppe (17.12%), and sandy vegetation (3.85%), and two other land use types: 
water body (0.57%) and Gobi desert (17.07%). 


3 Methodologies 
3.1 SPI 


In this study, we chose two drought indices to investigate GPP response to drought within the 
context of global warming. SPI, posited by McKee et al. (1993), only considers the influence of 
precipitation while disregarding the impact of temperature on drought. In 2010, Vicente-Serrano 
et al. (2010b) proposed the concept of SPEI, which accounts for the evaporation of precipitation 
due to temperature and thus acknowledges the influence of temperature change on drought. It offers 
a wider monitoring range for drought compared to SPI. Utilizing these two drought indices not 
only helps us distinguish the warming aridification and non-warming aridification areas but also 
quantifies the impact of warming on vegetation's response to droughts. 

SPI was calculated using monthly precipitation as the input data. In this study, the gamma 
distribution probability of precipitation was calculated and then normalized to find the output value 
of SPI (Neda et al., 2022). Assuming that the precipitation during a certain period is x (higher than 
zero), the probability density function of the gamma distribution is defined as follows: 


g(x)= xt B, (1) 


1 
PTa) 
T(a)= [, x7 tea , (2) 


where a (a>0) and £ ($>0) are the shape and scale parameters, respectively; (a) is the gamma 
function; and g(x) is the function for calculating the probability density function of the gamma 
distribution of precipitation. 

a and p can be obtained through the maximum likelihood method as follows: 


a- itvl+44/3 


pat (4) 
a 


A = In()-+ in, (5) 
i=l 
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where A is an intermediate variable; n is the length of the time series (month); x is the average 
value of precipitation at a certain pixel (mm); x; is the precipitation data series from CRU (mm); 
and 7 is the time series. Therefore, the probability of the random variable x being less than the 
precipitation xo in a certain year can be obtained as follows: 


x 1 x @-l-x/f 
G(x)= | g(x)dx = ee dx , (6) 
f 0 B°T( a) J 0 
where G(x) is the cumulative probability for a given time scale. 
Suppose that the precipitation during a certain period equals to zero, then the equation would 
be as follows: 


A(x)=q+(—q)G(x), (7) 

where q is the probability when precipitation equals to zero; and H(x) is transformed into a 
standard normal variable, which is the value of SPI. 

We used SPEI base v.2.7 dataset in this study. SPEI was calculated using precipitation and 

temperature as input datasets. Potential evapotranspiration calculation requires temperature data, 

and the difference between precipitation and potential evapotranspiration values was used as SPEI 


output value. The related calculation process follows the recommendations of Vicente-Serrano et 
al. (2010a, b): 


pet =16F( 7) (8) 

1.514 

T 

d'=|—| , 9 
e ° 
m=6.75x10/ I? —7.71x10° J? +1.79 x10 + 0.492, (10) 
oles aD 

12/\ 30 
N= (24 arecos -ian pran| 0.4093sen (2x -1 roj 5 (12) 

T 


where PET is the potential evapotranspiration (mm); T is the monthly mean temperature (°C); J is 
the heat index, calculated as the sum of 12 monthly index values d'; m is the coefficient associated 
with J; F is the correction coefficient computed as a function of latitude and month; NDM is the 
number of days in a month (d); N is the maximum number of solar hours (h); ø is the latitude (°); J 
is the average Julian day of the month (d); and 7 is the ratio of circumference to diameter. 
D=PRE-PET, (13) 
kel 
Di =)\(PRE,_;—PET,_;), z2k, (14) 
i=0 
where D is the SPEI; PRE is the precipitation (mm); « is the time scale of the aggregation; and z is 
the calculation number. 


3.2 Sen's slope estimation and Mann-Kendall (M-K) test 


In this study, Sen's slope estimation and M-K test were used to estimate the changes in vegetation 
GPP trends and various drought indices. Sen's slope estimation, a nonparametric statistical 
method, offers high computational efficiency and is insensitive to measurement errors and 
outliers. It can also effectively reflect the degree to which trends in a certain index over a long 
duration. Sen's slope (K) is calculated as follows (Agarwal et al., 2021): 


ia 
aera (15) 


52 JOURNAL OF ARID LAND 2024 Vol. 16 No. 1 


K =median(Q), (16) 
where Q is the slope; and Y; and Y; are the values of variable Y at times of i’ and j’, respectively 
(1<i’<j’<k). The median function sorts Q sequence from the smallest value to the largest value. If 
it is an odd sequence, the median value is taken as K value; and if it is an even sequence, the 
average of the middle two values is taken as K value. 

M-K test can be used to analyze and predict long-term climate series and is often employed to 
test the significance of changing trends by calculating Sen's slope. M-K test does not require data 
to follow a normal distribution and remains unaffected by missing values and outliers. Assuming 
the time series is {X;} (i=l, 2, ---, k) and satisfies the condition i’<j’<k, the statistics can be 
calculated as follows (Zhang et al., 2012; Kocsis et al., 2020): 


k-l k 
S=), >) sgn(X}-X)), (17) 
i'=1 j'=i'+1 
l, (Xj - Xj) >0 
sgn(X j — Xj.) = 50, (Xj -— Xp) =0 , (18) 
-l, (Xp -X;)<0 
E S>0 
aj var(S) 
Z= 50, S=0 ; (19) 
semis, S<0 
aj var(S) 


where S is the statistic being tested; X; and X; are the time series data for variables; Z is the 

standard normal test statistic, for a given significance level P, and |Z|>Z\_p2 implies that this time 

series has a significant changing trend at that level; and var(S) is the variance of S, which can be 

calculated as follows: 

k(k -1)(2k +5) 
18 l 


var(S) = (20) 


3.3 Hurst index 

Hurst index, based on the analysis method of the rescaled range, is widely used for monitoring the 
persistence of long-term time series changes (Li et al., 2021). In this study, this method was 
employed to monitor future continuous changes in vegetation GPP of the Mongolian Plateau. 
First, we divided the vegetation GPP time series (GPP) (t=1, 2, ..., 37) into t subsequences 


GPPw (1, 2, ..., T). GPP mean time series (GPP(,) ) was then calculated as follows: 


Se iE 
GPP) a aa (7 =1, 2, sees 37). (21) 
t= 


Subsequently, the cumulative deviation of pixel-by-pixel (GPPo,») was calculated using the 
following equation: 


i ETS 
GPP; a = >, (GPP -GPP >), (IS t<r). (22) 

The range sequence R,, was defined as follows: 
Ray = max GPP, 2) — min GPP; ,) (T =1, 2, ..., 37), (23) 


where maxGPPo,» and minGPPo,» (g C/(m?-a)) are the maximum and minimum values of GPP 
time series, respectively. 
The standard deviation sequence SD: was calculated as follows: 
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ee ; 
SD, = 7 (GPR = GPR) (24) 
t= 
R 
(7) H 
=(cT) , 25 
SD. (cr) (25) 


where c is the constant value; and H is the Hurst index. 
Hurst index was obtained using the following equation: 


where a is the constant value. 

Hurst index has three forms: if Hurst index is equal to 0.5, then time series changes randomly; 
if Hurst index is greater than 0.5 and less than 1.0, time series has long-term memory, and the 
future changes would be consistent with past changes; and if Hurst index is greater than 0.0 and 
less than 0.5, time series would be characterized by anti-persistence; in other words, the future 
persistent changes would oppose the past change trends (Qu et al., 2020). 

According to Sen's slope estimation and Hurst index analysis, we divided the sustainability of 
GPP changes in the Mongolian Plateau into six scenarios: a non-continuous increase (Hurst 
index<0.5 and Sen's slope>0.00), a non-continuous decrease (Hurst index<0.5 and Sen's 
slope<0.00), a continuous non-significant increase (Hurst index>0.5 and Sen's slope>0.00; 
P>0.05), a continuous non-significant decrease (Hurst index>0.5 and Sen's slope<0.00; P>0.05), a 
continuous significant increase (Hurst index>0.5 and Sen's slope>0.00; P<0.05), and a continuous 
significant decrease (Hurst index>0.5 and Sen's slope<0.00; P<0.05). 


3.4 Division of warming aridification areas 


Furthermore, in order to identify the contribution of temperature rise to drought, we divided the 
Mongolian Plateau into warming aridification areas and non-warming aridification areas based on 
the dry and wet change trends of SPI and SPEI. The criteria employed in this study were as 
follows: the areas indicated as not significantly drying when SPI was used but significantly drying 
when SPEI was used were identified as areas with mild aridification due to climate warming; the 
areas indicated as wetting when SPI was used but not significantly drying when SPEI was used 
were identified as areas with moderate aridification; and the areas indicated as wetting when SPI 
was used and significantly drying when SPEI was used were identified as areas with severe 
aridification (Zhang et al., 2022a). 


3.5 Correlation analysis and regression analysis 


In this study, the influence of temperature and drought indices on GPP was judged by correlation 
analysis and multiple correlation analysis. The correlation coefficient can represent the impact 
extent of climate factors on GPP, and the multiple correlation coefficient reflects the 
comprehensive influence of multiple independent variables on one dependent variable. 

Taking temperature as an example, the correlation coefficient (Rur) between GPP and 
temperature can be calculated using the following formula (Wu et al., 2022): 


Et Ur - DW, -V) 
Ruy = 7 = 7 = ; 
Era Ur -0 E0- 


where U and V are the average values of GPP (g C/m°) and temperature (°C), respectively; 
and U; and V; are the values of GPP (g C/m’) and temperature (°C) at time i', respectively. The 
value range of Rur is [—1, 1]. The larger the absolute value of Ruy, the greater the influence of 
temperature on GPP. 

We used the multiple correlation coefficient to identify the degree of influence of drought 
indices (SPI and SPEI) and temperature on GPP. The multiple correlation coefficient (ruw) can be 
calculated as follows: 


(27) 


54 JOURNAL OF ARID LAND 2024 Vol. 16 No. 1 


amw = 41- (1-2 -F wa). (28) 
where ru is the correlation between GPP and drought indices; and ruw,» is the partial correlation 
coefficient between GPP and drought index (SPI or SPEI) after the temperature is fixed. 

The regression coefficients can be used to compare and analyze the sensitivity of vegetation 
GPP to different climatic factors. Since the dimensions of the indicators are not uniform, they 
must be normalized using the following equation before regression analysis (Wang et al., 2018): 
E-E 
S(E) ` 
where Z(E) is the standardized time series; Æ is the GPP, SPI, SPEI, or temperature series; E is 
the average value of each series; and S(E) is the standard deviation of the corresponding series. 
The regression equation was then established, as follows, for the standardized parameters: 


Z(E) = (29) 


Yi =4;Xj7 +B, Xi9 + Ej» (30) 
where y; is the normalized GPP series; xip is the normalized temperature series; x/; represents the 
normalized drought indices; <; is the residual error; and a; and b; are the regression coefficients of 
temperature and drought indices, respectively. The values of a; and b; reflect the sensitivity of 
GPP to temperature and drought indices. 

In this study, we used "f-test" to determine the significance of the regression coefficient to 
ascertain whether the effect of the independent variables on the dependent variable was significant. 
The ratio of the regression coefficient and the standard error of the regression coefficient were used 


as t-values. Generally, P<0.05 represents a significant effect on GPP (Xu et al., 2003). 
3.6 Calculating CED 


CED refers to the cumulative effect of drought on vegetation on a continuous time scale. The 
correlation coefficient is calculated for monthly GPP values and drought indices over 1.00—12.00 
months. The cumulative month, where the maximum correlation coefficient is located, represents 
the cumulative effect of drought on vegetation GPP (Wei et al., 2022) as seen in the following two 
equations: 


r, = corr(GPP, AI) (1.00 </ <12.00), (31) 
Tnax—cum = Max(7;) (1.00 < Z <12.00), (32) 


where 7; is the correlation coefficient between GPP and drought indices; / is the monthly 
cumulative drought indices for 1.00—-12.00 months; AJ is the drought indices; and rmax-cum is the 
maximum correlation between GPP and drought indices. 


4 Results 


4.1 Trends in GPP and drought indices changes in the Mongolian Plateau 


Between 1982 and 2018, the changing rate of vegetation GPP in the Mongolian Plateau ranged 
from —10.64 to 12.42 g C/(m?.a), and the spatial heterogeneity in these trends was significant (Fig. 
2a). Vegetation GPP increased in 49.47% of the plateau's total vegetation-covered area (18.37% of 
which showed a significant increase), concentrated in the southeast and northwest regions of the 
Mongolian Plateau. Meanwhile, the areas with significant and non-significant downward trends in 
GPP accounted for 10.77% and 39.76%, respectively, primarily distributed along northern 
Mongolia and northeastern and central Inner Mongolia (Fig. 2b). Regarding the different 
vegetation types, the most notable GPP increase was observed in the restoration of sandy 
vegetation, with 66.09% of areas indicating a significant upward trend. It was followed by 
cropland and shrub, with 54.41% and 45.51% of the land, respectively, exhibiting an upward 
trend in GPP. In contrast, meadow steppe and deciduous forest experienced considerable 
degradation, with over 80.00% of these areas indicating a downward trend in GPP (Fig. 2c). 
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Fig. 2 Trends in gross primary productivity (GPP) changes and their persistence in the Mongolian Plateau. (a), 
spatial distribution of GPP changing rate; (b), classification of the significance of GPP changing trends; (c), GPP 
changing trends in different vegetation types; (d), spatial distribution of Hurst index; (e), persistence of GPP 
changing trends; (f), persistence of GPP changing trends in different vegetation types. NV, non-vegetation area; 
SI, significant increase; SDe, significant decrease; NSI, no significant increase; NSDe, no significant decrease; 
NCI, non-continuous increase; NCD, non-continuous decrease; CSI, continued significant increase; CSD, 
continued significant decrease; CNSI, continued non-significant increase; CNSD, continued non-significant 
decrease; CF, coniferous forest; MS, meadow steppe; SH, shrub; TS, typical steppe; DF, deciduous forest; AS, 
alpine steppe; CR, cropland; DS, desert steppe; SV, sandy vegetation. 


Combined with the results of Hurst index analysis, we also examined the sustainability of GPP 
trends in the Mongolian Plateau. GPP changes in 84.50% of the Mongolian Plateau showed a 
continuous and stable trend (Hurst index>0.5), of which 38.97% displayed a continuous upward 
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trend in GPP, mainly in the northern and southeastern parts of the plateau. In contrast, a 
continuous downward trend was mainly observed in the central part of the Mongolian Plateau 
(45.53%), suggesting a more promising outlook for vegetation changes in the southeast, 
northwest, and northeast marginal areas of the Mongolian Plateau than in the middle part. Areas 
of unstable changes were only noted in 15.50% of the total vegetation-covered area (Fig. 2d and 
e). Significant differences were observed in the persisting trends of GPP changes among different 
vegetation types. GPP of shrub, cropland, and sandy vegetation increased continuously (these 
areas had a continuous upward trend, accounting for more than 60.00% of these areas), while GPP 
of meadow steppe and deciduous forest decreased significantly (showing a continuous downward 
trend, accounting for more than 70.00% of these areas; Fig. 2f). 

The changes in SPI and SPEI in the Mongolian Plateau from 1982 to 2018 displayed strong 
spatial heterogeneity (Fig. 3). Since SPEI considers the effects of both precipitation and 
temperature changes, the drought area monitored using SPEI is wider than that monitored using 
SPI. Based on the results of changes concerning SPI, the area of the Mongolian Plateau indicating 
a drying trend accounted for 30.12% (Fig. 3b), which was far lower than the dried area monitored 
using SPEI (92.27%; Fig. 3e). This indicates an undeniable drying trend in the plateau caused by 
rising temperatures over the past 37 years. The drought-stricken areas of the Mongolian Plateau 
that were monitored using SPI were primarily concentrated in the Altay Mountains in the west, 
the Da Hinggan Ling Mountains in the east, and the Kente Mountains in the central region. 
SPI-based monitoring revealed that 69.88% of the Mongolian Plateau displayed a wetting trend 
(with 5.82% passing the significance test), mainly located in the southwest region (Fig. 3a and b). 
However, about 92.27% of the Mongolian Plateau showed a drying trend as per SPEI. Of these 
areas, 36.10% were identified as significantly drying, mainly located in the central and southern 
parts of the plateau. Only 7.73% of the areas exhibited an insignificant wetting trend; these areas 
were sporadically distributed in the eastern and northwestern mountainous areas (Fig. 3d and e). 
Significant discrepancies were evident in the findings after comparing SPI and SPEI assessments 
across different vegetation areas. SPI results indicated that typical steppe, desert steppe, cropland, 
and sandy vegetation areas were mainly wetting, while deciduous forest areas were drying in a 
large area (87.50%), indicating reduced precipitation in this area (Fig. 3c). According to 
SPEI-referenced findings, all the different vegetation types were mainly drying, especially in 
desert steppe (73.17%; Fig. 3f). 

By further dividing the observed drying and wetting trends associated with different drought 
indices, 71.39% of the Mongolian Plateau was found to be affected by varying severity of drought 
caused due to the rise in temperatures (Fig. 4). Of this area, 35.30% constituted moderate 
warming aridification areas, mainly concentrated in the Hangay Mountains in northern Mongolia, 
in western Inner Mongolia, and toward the south of the Da Hinggan Ling Mountains. The mild 
warming aridification areas closely corresponded to those areas indicated as not significantly dry 
based on SPI. Southern Mongolia and the central and western parts of Inner Mongolia were 
categorized as severe warming aridification areas, accounting for 28.53% of the entire plateau 
(Fig. 4a). Notably, desert steppe and sandy vegetation exhibited aridification due to warming, 
affecting 92.68% and 90.24% of the areas, respectively. Also, the severe warming aridification 
area in desert steppe accounted for more than half of its total land area (54.15%). Moderate and 
severe warming aridification areas accounted for 65.85% and 24.39% of sandy vegetation, 
respectively. Additionally, the arid areas of cropland and typical steppe accounted for more than 
50.00% of their territories, indicating that temperature rise greatly increases the drought stress of 
these vegetation types. In contrast, deciduous forest displayed the least sensitivity to temperature 
rise, with arid areas accounting for only 12.50% of their total area (Fig. 4b). 


4.2 Correlation analysis of GPP with temperature and drought indices 


The spatial distribution of correlation between different drought indices (SPI and SPEI) and GPP 
in the study area was similar, with primarily positive correlations (Fig. 5). The positive 
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Fig. 3 Variation trends of drought indices in the Mongolian Plateau. (a), spatial distribution of Sen's slope of the 
standardized precipitation index (SPI); (b), classification of changing trends of SPI; (c), SPI changes in different 
vegetation types; (d), spatial distribution of Sen's slope of the standardized precipitation evapotranspiration index 
(SPEI); (e), classification of changing trends of SPEI; (f), SPEI changes in different vegetation types. SH, 
significant humidification; SDr, significant dryness; NSH, non-significant humidification; NSDr, non-significant 


dryness. 


correlations of SPI and SPEI with GPP accounted for 89.66% and 81.98% of the total 
vegetation-covered area, respectively. Furthermore, significant positive correlations constituted 
31.71% and 21.45%, respectively, mainly concentrated in the eastern Mongolian Plateau. These 
results suggest that the vegetation in these areas is affected by moisture conditions. The areas where 
GPP exhibited negative correlations with SPI and SPEI were mainly distributed in the northwest, 
the eastern mountainous areas, and parts of the southwest, accounting for 10.33% and 18.02% of the 
total vegetation-covered area, respectively. Most of these areas were not significant (Fig. Sa, b, d, 
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Fig. 5 Spatial distribution and classification of correlation coefficients between GPP and SPI (a and b), GPP and 
SPEI (d and e), and GPP and temperature (g and h), and the percentage of the classification of correlation 
coefficients in different vegetation types (c, f, and i). SPC, significantly positive correlation; SNC, significantly 
negative correlation; NSP, non-significantly positive correlation; NSN, non-significantly negative correlation. 


and e). Meanwhile, GPP in these areas was mainly positively correlated with temperature, possibly 
due to the relatively high altitude, abundant moisture conditions, and low temperature of these 
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areas. Furthermore, the carbon sequestration of vegetation is mainly limited by thermal conditions. 
Additionally, a negative correlation was observed between temperature and GPP in other areas of 
the Mongolian Plateau (58.08%; Fig. 5g and h), indicating that increased temperatures did not 
correspond to an increase in the carbon sequestration capacity of local vegetation. 

GPP and drought indices of different vegetation types exhibited predominantly positive 
correlations (56.00%—95.00%). The positive correlation between GPP and SPI of cropland and 
deciduous forest was higher, accounting for more than 94.00% of the areas (Fig. 5c). The area 
with a positive correlation between GPP and SPI was larger than the area with a positive 
correlation between GPP and SPEI for various vegetation types (Fig. 5f). The statistical results 
regarding the correlation between GPP and temperature showed that GPP trends in meadow 
steppe, typical steppe, and desert steppe were in contrast to the changing trend of temperature. 
Consequently, an increase in temperature would negatively impact the growth of such vegetation. 
However, 78.28% of sandy vegetation GPP was positively correlated with temperature, indicating 
that temperature rise promotes GPP synthesis in sandy vegetation (Fig. 5i). 

The multiple correlation distributions of GPP with temperature and the two drought indices 
were relatively consistent (Fig. 6). The multiple correlation coefficients of GPP with SPI, SPEI, 
and temperature ranged from 0.00 to 0.81. The areas significantly affected by temperature and 
drought as per SPI accounted for 27.65% of the total vegetation-covered area, primarily 
distributed in the central, eastern, and southeastern parts of the Mongolian Plateau (Fig. 6b). In 
contrast, the area significantly correlated with temperature and SPEI showed a decrease. The 
decrease in significantly related areas was mainly concentrated in the Da Hinggan Ling 
Mountains (Fig. 6e). Notably, the multiple correlations among shrub, cropland, and sandy 
vegetation were significant, indicating that the combined action of temperature and drought has a 
stronger impact on GPP of these vegetation types (Fig. 6c and f). 


4.3 Sensitivity of GPP to drought and temperature 


An analysis of the sensitivity of GPP to drought indices and temperature from 1982 to 2018 
revealed that GPP in the Mongolian Plateau was more sensitive to drought than to temperature, 
especially in the central and eastern parts (Fig. 7). The regression coefficients of GPP with SPI, 
SPEI, and temperature ranged from —0.85 to 1.00 (Fig. 7a, b, e, and f). Overall, 65.20% (63.12%) 
of vegetation GPP in the Mongolian Plateau was found to be primarily and positively controlled 
by SPI (SPEI). Based on SPI and SPEI, the areas where GPP was primarily influenced by 
temperature accounted for 34.80% and 36.88% of the vegetation-covered area, respectively. The 
areas where temperature emerged as the most important limiting factor were similar, with SPI and 
SPEI as control, and were concentrated in the northwestern mountainous region and the 
southeastern region (Fig. 7c and g). The heightened positive sensitivity of GPP to drought, with 
SPEI and temperature as independent variables, was more pronounced than the correlation 
observed with SPI as independent variables (Fig. 7a and b). Since only the effect of precipitation 
is considered when SPI is the independent variable, while SPEI considers not only the effect of 
precipitation but also the effect of temperature, this suggests that the effect of temperature on 
drought increases the response of vegetation to drought sensitivity. Furthermore, the sensitivity of 
GPP to temperature and drought exhibited notable differences among various vegetation types. 
GPP of meadow steppe, typical steppe, deciduous forest, and desert steppe was found to be more 
sensitive to drought than to temperature, and the areas where drought was the primary control 
accounted for more than 60.00% (Fig. 7d and h). 

To determine the impact of warming on vegetation's response to drought, we calculated the 
sensitivity differences of GPP to drought in warming aridification areas and non-warming 
aridification areas across various vegetation types. The results indicate an alteration in the 
sensitivity of GPP to drought in various vegetation types due to rising temperatures (Fig. 8). 
Except for deciduous forest and sandy vegetation, the response of GPP to SPI increased with 
rising temperatures. Meanwhile, the response of coniferous forest, meadow steppe, deciduous 
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Fig. 6 Spatial distribution and classification of multiple correlation coefficients of GPP with SPI and 
temperature (a and b), and GPP with SPEI and temperature (d and e), and the percentage of the classification of 
multiple correlation coefficients in different vegetation types (c and f). SC, significant correlation; IC, 
insignificant correlation. 


forest, and alpine steppe to SPEI was aggravated, and the sensitivity of typical steppe, cropland, 
desert steppe, and sandy vegetation to SPEI was weakened. The difference in the sensitivity of 
GPP to drought between warming aridification areas and non-warming aridification areas of 
sandy vegetation was —0.24, which was considerably greater than that of other vegetation types. 


4.4 CED of drought on vegetation GPP 
Regarding SPI, the average CED of drought on vegetation GPP was about 5.77 months, with 


durations exceeding 10.00 months in the northeastern, southern, and central parts of the 
Mongolian Plateau (Fig. 9a). The average CED according to SPEI was 5.90 months. Overall, 
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Fig. 7 Sensitivity of GPP to drought indices and temperature and spatial distribution of main control factors. (a), 
spatial distribution of the sensitivity of GPP to SPI; (b), spatial distribution of the sensitivity of GPP to 
temperature with SPI and temperature as influencing variables; (c), spatial distribution of main control factors to 
GPP with SPI and temperature as influencing variables; (d), the percentage of main control factors to GPP in 
different vegetation types with SPI and temperature as influencing variables; (e), spatial distribution of the 
sensitivity of GPP to SPEI; (f), spatial distribution of the sensitivity of GPP to temperature with SPEI and 
temperature as influencing variables; (g), spatial distribution of main control factors to GPP with SPEI and 
temperature as influencing variables; (h), the percentage of main control factors to GPP in different vegetation 
types with SPEI and temperature as influencing variables. 
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Fig. 8 Effects of temperature rise on the sensitivity of GPP to drought in different vegetation types. (a), SPI; (b), 
SPEI. 
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CED of SPEI on vegetation GPP in the eastern region was longer than that in the northwestern 
region, indicating greater drought resilience in the eastern region (Fig. 9c). Based on SPI, the 
average CED of different vegetation types ranged from 5.02 to 7.33 months, with CED of 
deciduous forest and alpine steppe reaching 7.33 and 6.49 months, respectively; CED of other 
vegetation types remained concentrated around 5.00-6.00 months (Fig. 9b). Based on SPEI, the 
average CED for different vegetation types ranged from 3.03 to 8.95 months, with significant 
among different vegetation types. Coniferous forest, meadow steppe, typical steppe, alpine 
steppe, and desert steppe displayed relatively long CED, ranging from 6.78 to 8.95 months. In 
contrast, the drought resistance of shrub, deciduous forest, cropland, and sandy vegetation in the 
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Mongolian Plateau was relatively similar, with CED lasting for 3.00-4.00 months, which was 
comparable to SPEI (Fig. 9d). These observations imply that CED is closely associated with the 
differences in various vegetation conditions and the external environment. 


5 Discussion 


5.1 Sensitivity of vegetation GPP to drought 


The analysis of the temporal variation characteristics of vegetation GPP in the Mongolian Plateau 
from 1982 to 2018 reveals that the vegetation in the northwest and southeast regions of the 
plateau has recovered well over these 37 years, corroborating with previous findings (Meng et al., 
2023). This recovery shows the effectiveness and significance of the implemented ecological 
projects, such as sandstorm control and farmland-to-forest restoration. The results show that the 
drought severity in the Mongolian Plateau has significantly increased over the past 37 years as a 
result of rising temperatures, especially in typical steppe, cropland, desert steppe, and sandy 
vegetation areas. Considering that transpiration and evaporation consume around 80.00% of 
precipitation (Abramopoulos et al., 1988), SPEI is more suitable for monitoring changes in the 
drought severity of different regions, especially against the backdrop of global warming. 

The impact of climate warming on the drought sensitivity of GPP is more pronounced in 
coniferous forest, deciduous forest, meadow steppe, and alpine steppe. This is because 
vegetation growth in arid and semi-arid areas is mainly limited by insufficient moisture 
conditions or low water availability. This implies that, amidst global warming, changes in 
drought severity may have a greater impact on the role of vegetation as a carbon sink in these 
areas, which is consistent with the findings of extant research (Bai and Li, 2022). Rising 
temperatures lead to increased water loss via evapotranspiration, aggravating the severity of 
regional drought and resulting in more obvious water stress on vegetation (Li et al., 2018). The 
increase in the sensitivity of GPP to drought implies an increase in the intensity of drought 
effects on vegetation. Positive correlation was enhanced in areas with reduced water, while 
negative correlation was enhanced in areas with increased water, especially in the high latitudes 
of the northern hemisphere (Wei et al., 2023). 

It is worth noting that the temperature increase did not seem to augment the drought 
sensitivity of GPP in the deciduous forest and sandy vegetation when SPI was used for 
monitoring. However, the area primarily controlled by precipitation in the deciduous forest 
accounted for the largest proportion (78.27%), showing that this area is affected by changes in 
water conditions. Therefore, the changes in drought severity caused by temperature rise have not 
been obvious. However, the moisture conditions in sandy vegetation areas have improved 
significantly (Fig. 3), and the adverse effects of rising temperature on the vegetation in this area 
were slowed down by improved moisture conditions. Wang et al. (2022c) noted that in the past 
30 years, climate warming has led to warming and wetting trends in various areas of northern 
China, promoting the restoration of vegetation. Typical steppe, cropland, desert steppe, and 
sandy vegetation are mainly considered warming aridification areas (Fig. 4), and the sensitivity 
of their GPP to drought, as defined using SPEI, does not increase with climate warming (Fig. 8). 
Cropland and sandy vegetation GPP was more affected by temperature (Fig. 7), and only 37.68% 
and 24.39% of the vegetation areas were significantly dry, respectively. In high latitudes and 
humid areas in the northern hemisphere, the key factor limiting vegetation productivity is 
temperature rather than water availability (Ritter et al., 2020). So temperature warming did not 
increase the sensitivity of vegetation to drought in such areas. The vegetation productivity of 
typical steppe and desert steppe was mainly affected by precipitation, and the significant drying 
area of desert steppe reached 73.17%. This phenomenon may be due to the fact that vegetation 
productivity in being not only controlled by climatic factors is also influenced by vegetation 
physiology. For example, photosynthesis in the Amazon rainforest increases during the dry 
season (Green et al., 2020). 
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5.2 CED of drought on vegetation GPP 


In arid and semi-arid regions, moisture condition is the main factor controlling vegetation growth 
and development. The lagged impact of precipitation on vegetation significantly outweighs that of 
temperature, causing precipitation to have a more lasting impact (Nie et al., 2021). In this study, 
shrub, deciduous forest, cropland, and sandy vegetation were found to be more susceptible to 
short-term drought accumulation effects under SPEI metric. The main drought-causing factors 
vary based on the duration of the phenomenon. Short-term drought is mainly caused by a decrease 
in soil moisture content and river discharge in the source area, and medium-term drought is 
closely related to reservoir water storage capacity and river discharge in the middle reaches. 
Meanwhile, changes in groundwater reserves are a key factor causing long-term drought 
(Vicente-Serrano et al., 2010b). 

When compared with other vegetation types, the different types of grasslands in the Mongolian 
Plateau have longer average CED. This finding is consistent with previous studies indicating that 
grasslands in arid areas can adapt to long-term and high-intensity drought (Deng et al., 2020). 
Additionally, the difference in vegetation's root depth and biomass mostly explains the difference 
in CED. For example, grasslands primarily consist of herbaceous vegetation with shallow root 
systems, while forests comprise woody vegetation with a deep root system, allowing them to 
easily access water to maintain growth. Since it is generally known that CED of grassland is 
shorter than that of forest, the carbon sink of forest is more stable during droughts (Wolf et al., 
2013; Sun et al., 2015). However, this study's results show that the cumulative effect of drought 
on grasslands, as characterized by SPEI, is longer than in the case of deciduous forest, which 
means that the drought tolerance of the plateau's grassland is higher than its deciduous forest. This 
may be because these herbaceous plants grow in arid and semi-arid dry climates for a long time 
and are exposed to prolonged soil water stress. Vegetation can reduce water loss and improve 
water absorption capacity by changing its growth mechanism, thus enhancing its tolerance to 
long-term drought (Reddy et al., 2004). Additionally, the sensitivity of grasslands located on the 
desert periphery to changes in drought severity may not be higher than expected (Zhang et al., 
2022c). 


5.3 Vegetation GPP in response to drought at different altitudes 


Generally, as altitude increases, vegetation growth is influenced by a twofold dynamic: decreasing 
temperatures limit vegetation productivity while increasing precipitation promotes it (Piao et al., 
2006; Du et al., 2019). Therefore, this study further explored the influence of drought on 
vegetation GPP in the Mongolian Plateau at different altitudes. It is evident that at an altitude 
lower than 3500 m, the sensitivity of GPP to drought, when monitored using SPI and SPEI, 
decreased with an increase in altitude (Fig. Sla and b). These results highlight that the drought 
sensitivity of GPP is closely related to the long-term climatic conditions associated with regional 
vertical zonality. CED related to SPI decreased at first and then increased with altitude. At 
500-1000 m, the average CED was found to be the shortest (5.24 months), while the longest CED 
of 8.07 months was observed at an altitude higher than 3500 m. Overall, the cumulative effect of 
drought on GPP, as measured by SPEI, showed an upward trend, tending to stabilize at altitudes 
higher than 2000 m (Fig. Slc and d). This trend aligns with the findings of Zhang et al. (2022c), 
who noted that the cumulative effect of global drought and the drought tolerance of vegetation 
improved with an increase in altitude. 

Previous studies have reported that in 31.81% of drought-stricken vegetation-covered areas of 
the world, vegetation productivity is mainly affected by short-term drought (1.00-4.00 months) 
(Wu and Wang, 2022). The influence of drought accumulation time on vegetation in this study 
area was found to last much longer than the global average drought accumulation time, an 
observation that can be attributed to the longer CED of grassland at middle and high latitudes. For 
example, after an extended cold and dry period in northern Mongolia, which typically covers the 
land in ice and snow, the soil moisture in spring is mainly derived from the melting snow and 
thawing frozen soil, promoting the growth of grasslands (Nandintsetseg and Shinoda, 2011). It is 
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worth noting that different GPP estimation models may yield different vegetation photosynthesis 
simulations due to differences in their theoretical bases (Zhang et al., 2017; Wang et al., 2021a). To 
understand the response relationship between vegetation and drought in a more comprehensive and 
in-depth manner, it is necessary to compare and explore the phenomenon from multiple 
perspectives. 


5.4 Limitations and implications for future research 


A limitation of this study is that the Mongolian Plateau is dominated by agriculture and animal 
husbandry economy. Consequently, GPP is not only affected by climatic factors but also by 
human activities and other environmental factors, especially the persistent issue of overgrazing in 
Mongolia and central Inner Mongolia (Guo et al., 2021). Therefore, these factors should be taken 
into account in future studies. The data compiled in this study simulate the inter-annual variation 
of GPP and drought indices on a one-month scale, so studies of differences in vegetation GPP and 
dry and wet conditions in different months scales of a year could be added in future work. 
Additionally, due to limitations in the model's structure and parameter setting, this study was 
unable to consider a simulation of the relationship between the vegetation's physiological process 
and the environment; hence, the estimation accuracy of GPP could be further improved. Few 
studies in the extant literature focus on the variation of drought sensitivity with altitude, 
particularly in regions where altitudes surpass 3500 m. Therefore, future research can expand on 
the description on the relationship between drought and altitude presented in this study. 


6 Conclusions 


This study analyzed the spatiotemporal characteristics of vegetation GPP changes and examined 
their correlation with different drought indices in the Mongolian Plateau from 1982 to 2018, 
employing the long-term EC-LUE GPP and SPI and SPEI datasets. Overall, approximately 
49.47% of the total vegetation-covered areas within the region showed an upward trend in GPP, 
mainly concentrated in northwestern and southeastern Mongolia. The future GPP changes 
in 84.50% area of this region appeared to be continuous and stable. Among them, the 
restoration of sandy vegetation can be said to be the most significant, contrasted by the continued 
degradation of meadow steppe and deciduous forest. Warming aridification areas accounted for 
71.40% of the Mongolian Plateau. Severe warming aridification areas accounted for 28.53% 
of the region, mainly in the central and southern areas of the plateau. Around 54.15% of 
the desert steppe belonged to the warming aridification areas. The study confirms that 
climate warming has increased the drought sensitivity of GPP in coniferous forest, deciduous 
forest, meadow steppe, and alpine steppe. In addition, the plateau's vegetation 
demonstrated high adaptability to long-term drought, and the average CED for the whole 
region was around 6.00 months. It is noteworthy that under long-term water stress, CED, as 
monitored using SPEI, was generally higher for grassland than for deciduous forest. The 
drought sensitivity of GPP decreased with an increase in altitude. The average CED of GPP to 
SPI first decreased and then increased with altitude. The average CED per SPEI showed an 
upward trend at first and then stabilized. 
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Fig. S1 Altitude dependence of gross primary productivity (GPP) on drought sensitivity and drought 
accumulation. SPI, standardized precipitation index; SPEI, standardized precipitation evapotranspiration index. 
The boxes represent the range from the lower quantile (Q25) to the upper quantile (Q75). The cross symbols and 
horizontal lines inside the boxes represent the means and medians, respectively. The upper and lower whiskers 
indicate the maximum and minimum values, respectively. 


